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ABSTRACT 



This is a report of a seareh for propagation anomalies using a 
large quantity of high frequeney data produced as a byproduct of 
BRIGHAM, a Department of Defense projeel. The BRIGHAM data is 
based on 890 KHz wide samples of the HF spectrum at a 25 eycle rate, 
using a 2. 8 KHz resolution for a duration of approximately 2. 5 minutes. 
This method of data eolleetion is unique and it was hoped that propaga- 
tion anomalies, ineluding wide band anomalies, might be deteeted. 
Anomalies are believed to occur in the propagation of radio signals 
and they are usually other than known, routine effeets but may include 
known effeets which cannot adequately be explained. The seope of this 
examination was limited to the visual analysis of computer processed 
data presented on an interactive graphies unit. 



2 



TABLE OF CONTENTS 



I. OBJECTIVES 4 

II. INTRODUCTION 5 

A. NATURE OF THE DATA 7 

B. PRESENTATION OF THE DATA 7 

III. ANALYSIS 11 

A. PRELIMINARY 11 

B. PROCEDURE 12 

C. OBSERVATIONS 13 

IV. CONCLUSIONS AND RECOMMENDATIONS 36 

GLOSSARY 39 

COMPUTER PROGRAMS 40 

BIBLIOGRAPHY- 81 

INITIAL DISTRIBUTION LIST 2 82 

FORM DD 1473 83 



3 



I. OBJECTIVES 



This project provided the opportunity to examine unique data 
gathered with the AN/ELR-15 wide-band scanning receiver and the 
AN/ERD-10 circularly disposed antenna system. The project was 
suggested by the availability of this large volume of data which had 
been gathered at great expense for another purpose. The data seemed 
to offer great potential and the possible basis for considerable research 
and understanding about electromagnetic propagation. The initial con- 
cept was to observe data in the different form of a 3 -dimensional con- 
tour of amplitude, frequency and time. This concept included the use 
of a computer and interactive graphics unit to vary the size and dim- 
ensions of the contour to obtain a new view of the high frequency region 
of the spectrum and perhaps observe propagation characteristics not 
easily recognized otherwise. 

The objectives of this project were to take the given sets of 
data, write the computer programs necessary to present the data 
visually and to analyze the data. The objectives of the analysis were 
twofold: To study propagation effects on high frequency signals, 

becoming familiar with the data, the mode of visual presentation 
and the appearance of routine propagation effects; To look for unusual 



and anomalous propageition effects. 



II. INTRODUCTION 



It is stressed that the scope of this effort was confined to a 
visual analysis of the data for detecting unusual patterns and anomalies 
which could be attributed to propagation conditions. The search 
was conducted in the high frequency region (3-30 MHz) using data 
on computer cards. The data cards were produced during a phase 
of BRIGHAM, a project conducted for the Department of Defense by 
Sanders Associates, a private corporation. Each data set represents 
the output of a swept-tuned panoramic receiver and preserves amp- 
litude information vs frequency and time. For visual presentation, 
the data were processed on a Scientific Data Systems XDS 9300 
computer and displayed on an Adage AGT/10 graphics unit (figure 1). 
The equipment is part of the Electrical Engineering Computer Center 
of the Naval Postgraduate School. Interesting sets of sweeps have 
been extracted as one picture, plotted on a CALCOMP 563 plotter 
and included in this report. Forty-six data sets were available for 
examination, each consisting of 7195 cards. Several data sets were 
examined in detail with the intent of using other sets for verification 
of observed anomalies. Recommendations for further examination 



are given at the end of this report. 







FIGURE I s THE VARIOUS STAGES OF PROCESSING DATA 
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A. 



NATURE OF THE DATA 



The data were gathered using a swept-tuned receiver and either 
an omnidirectional or a 12 beam antenna. Each data set was produced 
by sampling the receiver output as it scanned downward through an 
890 KHz band. Sets were gathered at different base (starting) fre- 
quencies and at different times of the day. The sets contain the 
amplitude (in db) of the 318 contiguous frequency bins (sampling points) 
and consist of 7194 cards to represent the 3597 scans. The initial 
card in each set contains date, time and parameter information (fre- 
quency, etc. ) for that set. 

B. PRESENTATION OF THE DATA 

Computer programs for (1) transferring the data from cards to 
tape, (2) visual presentation of the data and (3) extraction of the data 
for graph plotting are included at the end of this report. There are 
two computer programs for presenting the data and they reflect the 
trade-off between the width of the displayed spectrum and the number 
of scans which could be viewed simultaneously. This restriction was 
due to core memory size in the XDS 9300 and again in the AGT/10. 

1 . Anomaly A 

This program was used for viewing the data as a contour 
and required appreciable dimensions in time and frequency to be of 
use. It was initially adapted from a library program but, along with 
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two accompanying subroutines, was changed almost entirely in adapt- 
ing it to the BRIGHAM data. Anomaly Plot A was used with this 
program for plotting graphs. The data were displayed as shown in 
figure 2 with the most recent scan at the bottom of the picture. For 
each picture change, a new scan was brought on at the bottom and the 
oldest scan moved off at the top. The scans were not moved up in 
a continuous or movie-like fashion but all shifted upward a step at a 
time to replace previous scans. Options were available for bringing 
on more than one scan at a time. The amount could be varied so as 
to change the entire picture at once if desired. Scans could be skipped, 
taking every fifth or tenth scan for example, effectively decreasing 
the receiver sampling rate to accentuate long term or slowly changing 
trends in the data. Also, the time axis could be rotated to vary the 
3-dimension or surface contour effect. Other features are described 
in the computer program section. An optimum size of 130 bins and 
20 scans was chosen for this program. A deficiency in this presenta- 
tion was that amplitude changes with time were not easily observed. 

To overcome this, the Anomaly B program was written. The scale 
at the bottom of figure 3 and all the graphs is 10 frequency bins 
(28 KHz) per division. The leftmost bin in each group is labeled with 
the computer variable ISTRT, for the integer value of the "starting" 
bin. ISTRT is a NAMELIST variable and its value can be changed 
at the graphics console. It enabled the viewer to look at different 
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portions of the band. For identification, bins are located by counting 
to the right from the starting bin (frequency decreasing toward the 
right). The vertical scale of the Anomaly A presentation is 5 db of 
amplitude per division (referred to receiver sensitivity). Scans are 
spaced 1/25 second apart in time or multiples thereof when scans are 
skipped. 

2. Anomaly B 

This program was used to examine more carefully the 
changes in amplitude among a number of signals (see figure 3). At 
the top are five traces of signals chosen by specifying the bin number. 
The traces show the amplitude history of each signal up to a maximum 
strength of 30 db each, beginning with the present base scan at the 
left and extending to the right, showing the last 80 scans. The am- 
plitude was displayed after taking the maximum value of the signal 
in three bins (the one selected and one on either side) and scaling 
it. At the bottom are displayed all 318 bins to aid in selecting the 
five signals to be viewed. Because of the memory size trade-off, 
only five scans are shown simultaneously. The vertical scales 
differ from those in Anomaly A. Because of a smaller scale, the 
divisions are 10 db apart. Anomaly Plot B was used to plot this 
display on paper. 
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III. ANALYSIS 



A. PRELIMINARY OBSERVATIONS 

Before beginning the visual analysis, it was useful to predict 
the type of displays the data could yield. These preliminary observa- 
tions helped visualize numerous features and signal variations not 
due to propagation so that they could be recognized and not be inter- 
preted as anomalies to be checked and verified. Examples of these 
features included the following: 

1. Transmitter characteristics: Turning on and off; different 

types of keying which affect frequency and/or amplitude. 

2. Receiver characteristics: Filters for notching out inter- 

fering signals. 

Other preliminary observations predicted limitations on the 
ability to observe valid propagation characteristics which might be 
present. These limitations were due to receiver characteristics, 
the nature of the data and the method of visual presentation. 

1. Limitations of the receiver: The 890 KHz band restricted 

the ability to relate the effects on harmonics in the megahertz region 
and the ability to observe variations over a wide band. 

2. Limitations due to the data: The lack of resolution 

imposed by the 2. 8 KHz wide frequency bins was expected to make 
it impossible to measure slight variations in frequency. Better 
resolution could provide better identification of the types of signals 
present for the purpose of selecting ones with more stable transmitted 
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characteristics to permit more reliable analysis and measurement. 
The short duration (2. 5 minutes) sample provided a mountain of data 
but was expected to prevent observation of some propagation charac- 
teristics which require a longer period to manifest themselves. 

3. Shortcomings in the visual mode of presentation were not 
foreseen until much later in the project. At first only a general con- 
cept of the type of presentation existed and the final characteristics 
were considerably different. In fact, Anomaly B was developed 
only after much analysis and frustration with Anomaly A. It should 
also be noted that although the B version was developed last, it was 
not a replacement and both programs complement each other. 

B. PROCEDURE 

After reading a data set onto tape, the data were analyzed 
using the Anomaly A presentation to examine the entire length and 
width of the set. The object at this stage was to find gross anomalies 
in the 3 -dimensional contour which might be easily recognized. 
Hypothetical examples hoped for might have been any of the following: 
A rapid change in one region only of the contour; similar changes to 
one or more amplitude humps (signals or combinations of signals) in 
separated but perhaps related areas of the contour; holes or the 
absence of signals in limited bands of the contour. The Anomaly B 
program was then used for more c ritical analysis of amplitudes over 
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a longer period* The amplitude-time traces displayed data over a 
period four times as long as in Anomaly A (80 data points vs the 20 
scans) and examined it from a different perspective (abscissa in 
dimensions of time rather than frequency to examine individual 
signals ). 

C. FINDINGS AND OBSERVATIONS 

Data sets 4, 5, 8, 9, 15, 32, 34 and 147 were examined in 
detail. Set 4 was the first set examined and received considerable 
attention while being used to develop the first computer program. A 
routine for numbering each scan and later displaying the number of 
the base or bottom scan on the screen was developed and became 
indispensable in helping identify the location of anomalies or unusual 
data for subsequent observations. 

The first unusual feature noted is shown in figure 4 (upper left 
hand corner). The amplitude hump there is about 4-6 KHz wide and 
is seen to move -downward in frequency with time at a rate of about 
250 KHz per second. An effort to trace it through other signals was 
unsuccessful. The left hand margin represents the upper limit of 
the band (18 MHz), therefore tracing it at higher frequencies was not 
possible. The rate of frequency change was too fast to be classified 
as a ''whistler 11 or "dawn chorus". These phenomena are only found 
at considerably lower frequencies (several kilohertz) and the bandwidth 
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FIGURE 4. Data set 4 
ISTRT 1 
Base scan 2109 
Base frequency 18 MHz 
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of the amplitude hump was too great for this type of classification. 

This was the only occurrence of this phenomenon which was observed. 
It may have been a transmitter being tuned or a special purpose signal. 

A feature present in sets 4, 5, 8, 9, 15, 34 and 147 are sudden 
discontinuities across the entire spectrum of the set. Examples are 
shown in figures 5 through 7. The transitions were too rapid to be 
attributed to propagation effects and a satisfactory explanation has 
not been found. Late in the analysis, in set 34, several pictures from 
Anomaly A were pieced together (figure 7). The figure shows an 
effective shift of 160 scans. At one point it seemed to be the result 
of a missing data card. Two cards are required to describe one scan 
and if one is missing, the next card and all subsequent cards are all 
shifted back 160 scans. It would give the appearance on the left half 
of the scan being exchanged with the right half, as it appears in 
figure 7. A second missing card would then put them back in the 
right order, ending the "discontinuity”. Unfortunately, this simple 
explanation was not correct because an examination of the data in the 
region of the set 34 discontinuity revealed no missing or out-of-place 
cards. Examination of cards in other sets with discontinuities also 
revealed no missing cards. The bad data therefore must have been 
due to some sort of receiver or recording malfunction. 

Figure 8, instead of showing any anomalies, shows a set of 
well ordered, continuous signals of nearly constant amplitude and 
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FIGURE 5. Data set 4 
ISTRT 1 
Base scan 1423 
Base frequency 18 MHz 
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FIGURE G. Data set 4 
ISTRT 100 
Base scan 949 
Base frequency 18 MHz 
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e frequency 5 MHz 
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shape. This is contrasted in figure 9 by an irregular contour. In 
addition to variations in amplitude, some signals seem to be turned 
on and off alternately. The region between bins 30 and 38 contains 
a signal regularly occurring only once in nearly every three scans. 

At bin 60, the pattern is irregular and is probably an on-off keyed 
signal. 

When figure 10 was first observed, it was thought that a genuine 
anomaly has finally been found. Band limited noise seemed to appear 
first to the right of bin 160 and then spread to higher frequencies 
before subsiding. The presence of the notch seemed connected also, 
marking the lower limit during the first scan and persisting before, 
throughout and after the noise in time. After some study and con- 
sideration it was realized that each sweep required a finite amount 
of time and that the onset of noise appeared, although abruptly, at 
midscan. It is seen to grow weaker with time before finally 
disappearing. 

Further analysis was conducted with Anomaly B to examine 

amplitude variations of individual amplitude humps. Examination 

was begun in set 32 because it was free of the discontinuities present 

o 

in each of the previous sets. Also, set 32 was taken from a 12 beam 
antenna rather than the omnidirectional one used in the sets from 
4 to 15. The narrow beam was chosen to confine the signals observed 
to a particular region of origin. Although points of signal origin, 
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FIGURE 9. Data set 8 
ISTRT 1 
Base scan 60 
Base frequency 5 MHz 
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FIGURE 10. Data set 8 
ISTRT 100 
Base scan 1137 
Base fi'equency 5 MHz 
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and consequently propagation conditions, would still vary with each 
signal, the variations would not be as extreme as in the 360° case. 

Figures 11, 12 and 13 were extracted from set 32 and compare 
the amplitudes of signals found in bins 19, 64, 101, 156 and 182. 

Each of the three figures were drawn with a different sampling rate 
but have the same terminal or base scan (leftmost points at the top 
of the figures and the bottom scan in the lower section of the figures). 
Figure 11 shows every other scan (the variable NSKIP set equal to 1 
in the computer program) beginning with scam 2840, covering a 6. 4 
second period. Figure 12 shows every tenth scan (NSKIP^IO) begin- 
ning with scan 2200 (covers 32 seconds). Figure 13 shows every 35th 
s can beginning with scan 200 (covers 112 seconds). The scans at the 
bottom of each picture represent the entire data spectrum and show 
the last five scans at the particular skip rate noted. 

A few similarities exist between the traces in figures 11 through 
13. Propagation effects had been expected to create corresponding 
dips in amplitude among the signals of similar points of origin (such 
as those shown) but they were expected to appear over longer periods. 
On the other hand however, the rapidity of the fluctuations is not 
unusual. Ionospheric surface irregularities and layer displacements 
could account for them as described in reference 1. David and Voge 
[Ref. 2], for example, report that horizontal displacements have been 
measured by doppler techniques at velocities up to 60-80 miles per 
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FIGURE 12. Data set 32 
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hour. At the top right in figure 11, traces 2 and 3 show an inverse 
relationship. Traces in figure 12 indicate that the higher frequencies 
were affected first with the effects shifting lower with time. It is 
realized that since the signals selected for observation here were 
not equally spaced in frequency, the straight line analysis may not 
be particularly valid. Also, because of the small number of similar- 
ities in these samples, they may only be coincidental. 

Figures 14 through 17 were also taken from set 32 but compare 
the amplitudes of bins 35, 85, 167, 224 and 303. The sampling rates 
are shown for each figure as well as the terminal scan number and 
duration of the sample. Figure 14 shows no anomalies or correlation 
among signals with NSKIP=2. In figure 15, a periodic relation of longer 
duration seems to exist between traces 1 and 3. The propagation con- 
ditions affecting the two signals are similar but with the higher fre- 
quency signal (trace 1 at the top) being affected earlier in time. These 
effects are seen again in figure 16 but at a higher skip rate (dips 
closer together) where they can be observed over a longer period. 

No data exists to verify this pattern for the preceding time interval 
but it does not appear in the following interval (figure 17 - covers 
the interval from scan 1600 to number 3200 - NSKIP=20 over 80 data 
points). It may have died out or the display characteristics may have 
prevented it from being seen. The relationship stands out because 
of the frequency dependence of the effects. This type of anomaly is 
characteristic of what was expected when the project was begun. 
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Base scan 1000 trace 2 bin 85 

Base frequency 11 MHz trace 3 bin 167 

NSKIP 6 trace 4 bin 224 

Sample duration 19.2 sec trace 5 bin 303 
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FIGURE 16. Data set 32 
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FIGURE 17. 
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Figure 18 shows an unusual relationship between the signals of 
traces 1, 2 and 3. Two of the signals experienced some sort of- 
attenuation while the signal in the middle was enhanced. A similar 
inverse relationship was observed in figure 11. The same phenomenon 
as in figure 18 can be seen in figure 19 but at a higher skip rate. Also 
present were periodic variations in amplitude, particularly in trace 3. 
Periodic variations became more pronounced in trace 4 (figure 20) 
in a later time interval. 

In addition to observing those groups of signals already noted, 
other groups of signals were also chosen in order to sample a wide 
range of possible combinations. Some groups spread the sample over 
the entire data spectrum while others were confined to a narrow range 
of bins. A spread of low and high skip numbers was also included. 
Specifically, the following groups of signals in set 32 were among 
those examined but with nothing significant observed: Bins 3, 9, 16, 

21 and 25 with NSKIP equal to 1, 4, 16 and 40; Bins 135, 147, 158, 

168 and 177 with NSKIP equal to 3, 12 and 45; Bins 8, 13, 17, 22 and 
27 with NSKIP equal to 2, 10 and 35. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



Although a tremendous amount of data was available as a by- 
product of project BRIGHAM and a large quantity of it was examined in 
detail during this effort, the many anomalies 'hoped for at the begin- 
ning were not seen. Three significant phenomena were observed and 
preserved in this report. The first, the signal decreasing in freq- 
uency with time was almost certainly not a product of a propagation 
anomaly. Its identification might prove to be an interesting challenge, 
however. The second phenomenon showed that its periodic features 
were frequency dependent. Although periodic fading or intensity 
fluctuations were noticed at various times while examining the data, 
none were so well defined nor were they so well related between 
signals as in this case. Comparing the ionospheric layers to a cloud 
whose features are constantly changing, it was not surprising to see 
the phenomenon confined to a short period of time. Likewise, the 
irregular features of the ionosphere would allow several patterns 
to occur at once, superimposed on each other to create the overall 
observed effect. Some of these secondary features seem to show 
through in figure 16. The third phenomenon was the simultaneous 
enhancement and attenuation of signals. The seemingly incompatible 
conditions which caused this may have only been due to multipath 
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effects. It is strange, however, that no such pronounced effects 
were observed before or after the anomaly, particularly in trace 3, 
and that they occurred simultaneously among the three signals. It is 
expected that a great number of other anomalies exist and that many 
of them are preserved in the BRIGHAM data. Approximately one 
fifth of the available data were examined and due to time limitations, 
only one fifth of that amount was examined with Anomaly B. Analysis 
with Anomaly B was particularly time consuming because only five 
signals could be examined at once. It is likely that other anomalies 
are observable with the existing programs and many others detectable 
with modifications and different analysis tools. Some recommendations 
for improving the analysis techniques are as follows: 

1. Amend the Tape Load program to read the binary number 
of each data card and alert the operator to missing cards. A one-up 
numbering system to keep pace with card numbers could be used 
with an "if 11 statement to trigger an output statement when a card is 
discovered missing. 

2. Modify the procedure for selecting data for the amplitude 
history traces in Anomaly B. When cards are skipped to examine 
long-term propagation effects, an averaging routine should be entered 
to remove the short term amplitude fluctuations. The presence of 
the short term effects tends to mask those manifested only over a 
longer period of time. A modification of this type was attempted 
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near the end of the analysis but was unsuccessful due to a lack of 
time needed to perfect it. 

3. Incorporate a feature for backing up the tape by a variable 
(NAMELIST variable) number of scans. Two attempts to do this 
failed. In one case it was due to the fact that each scan constituted 
one '’physical record" in length and that backspacing was based on 
"logical records". If accomplished, the feature would be a significant 
time saving addition. 

4. Temporarily abandon the visual analysis approach of 
Anomaly B and employ statistical analysis techniques. Evaluation 
of the correlation between various sets of signals could be conducted 
rapidly and rather exhaustively on a larger computer such as the 
IBM 360 available at the Naval Postgraduate School. Preliminary 
analysis with Anomaly A might prove useful in the selection of signals 
for analysis and indispensable in locating the discontinuities described 
earlier so as to avoid misinterpreting results. 

5. Use of a larger, faster computer such as the IBM 360 to 
load cards onto magnetic tape would significantly reduce "dead time" 
in an extension of this project. Over an hour was required to load a 
single data set onto tape with the equipment described in this report. 
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GLOSSARY 



1. Bin (frequency bin): A 2. 8 KHz wide data sampling point. 

Energy present in these limits was quantized and represented in the 
range from 0-28 db. 

2. Data set: 3597 successive sampled sweeps or scans of the 

receiver and preserved on punched cards of magnetic tape. 

3. Dawn chorus or whistler: Interference in the audio range from 

storm discharges or high energy particle interaction in the ionosphere. 
Dawn chorus is produced around sunrise and sounds like a serenade 

of birds. It is a rise in frequency above about 2 KHz. Whistlers are 
audible decreases in frequency and are propagated from storm dis- 
charges along magnetic force lines of the earth in the range from a 
few hundred hertz to 20-30 KHz. 



4. NAMELIST variables: Computer variables, properly specified, 

which can be entered or changed after the program has been compiled 
and is running. These variables were used to control various aspects 
of the program. The most used were: 



ISTRT - 

IWDTH _ 
NSKIP - 



Integer number of the leftmost or starting bin to be 
displayed. 

Integer number of adjacent bins to be displayed. 
Number of scans to be skipped over. 



NMBRJMP - Number of scans to be brought on (and moved off) 
the screen at one time. 

5. Scan: One sweep of the receiver from one end of the band to the 

other. 
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140 CQNTINUE 
RETURN 
END 



SUBR9UTINE D I SPLAY ( NSCAN* I WDJH ) 

INTEGER AXES ( 47 ) 

COMMON /AREA1/FILE( 130,20 ), MOV (20/ 130), I MAGE (2700) 

COMMON /AREA2/IDEV/XSLNDYSLNT/X1/ YUDX*AXES#NMBRJMP/NRSCAN 
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The following subroutines, GINP, GINPUT, FNS and VCD, are 



common to and required by both programs Anomaly A and Anomaly B. 
They are library subroutines but are nonstandard and subject to change. 
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SU3R8UT I NE G I NP ( IDEV*ITDIR> I BLK/ IBUF) 
DIMENSION I BUF ( 1 ) * 1 TD I R ( 1 ) 



SUBROUTINE "GINP" 
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SUBROUTINE "GINPUT" 
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SUBROUTINE "FNS" 
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SUBROUTINE "VCD" 
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QUIRES SENSE SW« 2 AND r 0|_ C0NTR9L CARDS (LEFT JUSTIFIED T0 C0L 1) 
SPATCH 
$>>DATA 
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PROGRAM "GRAPH PLOT' 



"GRAPH PLOT M is a library program on paper tape and is used with 
the CDC 160 computer to drive the CALCOMP 563 plotter. It is 
found in the Electrical Engineering Computer Center of the U. S. 
Naval Postgraduate School. 
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